Nota: Esta es una aproximación cruda al modelo epidemiológico, y no debería tomarse como un resultado oficial ni de alta confiabilidad. Su interés es meramente referencial.
Usaremos la ecuación de Gompertz:
\[ y(t) = a \times e^{-e^{\frac{(\mu \times e)}{a} \times (\lambda - t) + 1}} \] Donde los parámetros a ajustar son a, mu (\(\mu\)) y lambda (\(\lambda\))
El resultado de usar todos los datos disponibles se ven a continuación.
Formula: y ~ a * exp(-exp(mu * exp(1)/a * (lambda - t) + 1))
Parameters:
Estimate Std. Error t value Pr(>|t|)
a 938.4392 226.4355 4.144 0.000992 ***
mu 49.6464 4.1250 12.036 9.0e-09 ***
lambda 8.8066 0.5512 15.979 2.2e-10 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 11.66 on 14 degrees of freedom
Number of iterations to convergence: 8
Achieved convergence tolerance: 6.266e-06Además, esta ecuación tiene un valor de AIC de 136.4641714.
El gráfico que acompaña a este modelo, presenta (en rojo) la extrapolación a 7 días de la curva actual. Además, se muestra, en color naranja, la curva extrapolada usando los datos de los primeros 12 días desde el primer caso confirmado (que fue el 20202-03-06).
Esta visualización esta hecha usando el lenguaje de programación R, y hace uso de una serie de librerías (flexdashboard, ggplot2, echarts4r, DT), para producir los distintos gráficos y tablas.